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ABSTRACT 

We have undertaken a survey of N 2 H + and N 2 D + towards 31 low-mass starless cores using 
the IRAM 30-m telescope. Our main objective has been to determine the abundance ratio of 
N 2 D + and N 2 H + towards the nuclei of these cores and thus to obtain estimates of the degree of 
deuterium enrichment, a symptom of advanced chemical evolution according to current models. 
We find that the N (N 2 D + )/7V (N 2 H + ) ratio is larger in more "centrally concentrated cores" with 
larger peak H 2 and N 2 H + column density than the sample mean. The deuterium enrichment in 
starless cores is presently ascribed to depletion of CO in the high density (> 3 x 10 4 cm~ 3 ) core 
nucleus. To substantiate this picture, we compare our results with observations in dust emission 
at 1.2 mm and in two transitions of C ls O. We find a good correlation between deuterium 
fractionation and iV(C 18 0)/iV(H 2 )i. 2 mm for the nuclei of 14 starless cores. We thus identified 
a set of properties that characterize the most evolved, or "pre-stellar" , starless cores. These 
are: higher N 2 H + and N 2 D + column densities, higher 7V(N 2 D + )/7V(N 2 H + ), more pronounced 
CO depletion, broader N 2 H + lines with infall asymmetry, higher central H 2 column densities and 
a more compact density profile than in the average core. We conclude that this combination 
of properties gives a reliable indication of the evolutionary state of the core. Seven cores in 
our sample (L1521F, OphD, L429, L694, L183, L1544 and TMC2) show the majority of these 
features and thus are believed to be closer to forming a protostar than are the other members of 
our sample. Finally, we note that the subsample of Taurus cores behaves more homogeneously 
than the total sample, an indication that the external environment could play an important role 
in the core evolution. 
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1. Introduction 

Understanding the physical and kinematic con- 
ditions in a molecular cloud core just prior to the 
onset of gravitational collapse is of paramount im- 
portance if one wishes to study the subsequent 
protostellar evolution. Molecular line measure- 
ments can provide the kinematical informations 
to probe infall motions in the highest density nu- 
cleus as long as the observed species is reasonably 
abundant in the region of interest. 
An important step forward in starless core stud- 
ies has been the understanding of the fact that, at 
temperatures of roughly 10 K and densities above 
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a few times 10 4 cm" 3 , CO, as well as other carbon- 
bearing molecules, condenses out onto dust grain 
surfaces preventing us from observing the high- 
density structures from where one expects the first 
signs of collapse (Leger 1983). This is most clearly 
seen from the comparison of maps made in the mil- 
limeter dust continuum with maps of CO isotopo- 
logues (Kuiper, Langer, & Velusamy 1996; Kramer 
et al. 1999; Alves, Lada, & Lada 1999; Casein et 
al. 1999). 

For reasons that are thought to have to do with 
the low binding energy of the N 2 molecule on ice 
mantle surfaces (Sadlej, Rowland, Devlin, & Buch 
1995; Bergin & Langer 1997), there are some N- 
containing molecular species that survive in the 
gas phase at least for densities in the range 10 5 - 
10 6 cm~ 3 . In particular, maps of N 2 H+ and NH 3 
(as well as their deuterated counterparts) show 
reasonable agreement with the distribution of dust 
emission. Thus, only from observations of such 
species one can hope to probe the physical and 
chemical conditions of the most evolved pre-stellar 
cores that show a high column density peak in the 
dust emission. 

One can also hope that by studying samples of 
such high density objects it will be possible to 
identify the cores where collapse has already been 
initiated or is at the point of doing so. 

The aim of the this paper is to report the results 
of a survey of 31 cores in various transitions of 
N 2 H+ and N 2 D+ observed with the IRAM 30-m 
telescope. 

These observations, besides probing the kine- 
matics of the high density nucleus of starless cores, 
enable us to gauge the deuterium enrichment in 
our sample. This parameter, evaluated from the 
N 2 D + over N 2 H + column density ratio, can be 
used as a tool to discriminate the evolutionary 
status of a starless core. In fact, the disappear- 
ance of CO from the gas phase, together with 
the exothermicity of the production of H 2 D + , 
leads to abundances of deuterated species in cold 
molecular clouds more than 3 orders of magni- 
tudes larger than the [D]/[H] elemental abundance 
(~ 1.5 x 1CT 5 Oliveira et al. 2003). The freeze- 
out of the CO onto the grain mantles decreases the 
destruction of H 2 D+ and . As a result, the en- 
hanced Hg abundance speeds up the reaction that 
leads to a higher production of H 2 D + and this in- 
creases the ratio of [H 2 D+]/[H+] (e.g. Dalgarno & 



Lepp 1984; Roberts & Millar 2000). 

We expect to use the N (N 2 D+)/7V(N 2 H+) ratio 
to discriminate the less evolved starless cores from 
the more advanced pre-stellar cores. For practi- 
cal purposes in this study, we define "pre-stellar 
cores" to be those with an A(N 2 D+)/7V (N 2 H+) 
abundance ratio > 0.1, and we believe that this 
approximately corresponds to a core with cen- 
tral density greater than 5 x 10 5 cm~ 3 that is in 
the last phase of its evolution prior to dynami- 
cal collapse. We base this partly on the cross- 
correlation of our iV(N 2 D + )/A(N 2 H+) measure- 
ments with other properties expected to probe an 
enhanced evolution of the starless cores such as 
broader line widths, larger n(H 2 ) central density, 
or larger amount of CO freeze-out. For this rea- 
son, supplemental dust continuum and C 18 mea- 
surements, able to furnish the H 2 density profile 
and the CO depletion, were carried out in cases in 
which these were not available in the literature. 

In § 2 of this paper, we outline our observational 
technique both for the N 2 H + and N 2 D + observa- 
tions and for the dust continuum and C 18 . In 
§ 3, we summarize our direct observational results, 
and in § 4 wc derive physical and chemical quan- 
tities for the cores in our sample. In § 5, we dis- 
cuss statistical correlations between the physical 
properties derived for our sample and in § 6, we 
summarize our conclusions. 

2. Observations 
2.1. The Sample 

We have selected 31 cores to be observed in 
N 2 H + and N 2 D + . The criteria used in the sample 
selection were not completely homogeneous, but 
we think that the sample should be representative 
of evolved starless cores. 

"Evolved" in this context means that there 
is evidence for high density in the form cither 
of strong N 2 H + (l-0) emission (> 3 K km s _1 
in Tmb scale) or strong 1.2-mm dust continuum 
(> 30 mJy/(ll" beam)) emission or both. The 
rationale here is that both of these criteria sug- 
gest a high central H 2 column density and hence 
high central pressure, given that N 2 H + appears to 
trace the dust well (e.g. Tafalla et al. 2002). 

Of the 31 objects selected in this manner, 24 
belong to an earlier survey of N 2 H+ and CS by 
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Core Sample 
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Note. — Col (1) most common core name; col. (2) and (3) reference position in J2000; col. (4) 
distance; col. (5) references for the position and distance, see below for coding; col. (6) Star 
forming environment: T: Taurus molecular cloud; O: Ophiucus molecular cloud; B: Bok globules; 
A: Aquila rift; col. (7) number of observed positions in both N2H+ and N2D+; col. (8) availability 
of a 1.2-mm continuum map and reference to it, coding below; col. (9) availability of a C ls O map 
and reference to it, coding below. 

References. — a: Lee, Myers, & Tafalla (2001); b: Tafalla et al. (2002); c: Caselli, Benson, 
Myers & Tafalla (2002); d: Bacmann et al. (2000); c: Shirley, Evans, Rawlings, & Gregersen 
(2000); f: Bergin, Alvcs, Huard, & Lada (2002); g: Crapsi et al. (2004); h: this work; i: Ward- 
Thompson, Motte & Andre (1999); 1: Pagani et al. (2003); m: Bianchi et al. (2003); n: Tafalla et 
al. (in preparation); o: Caselli et al. (1999); p: Bacmann et al. (2002). 
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Lcc, Myers, & Tafalla (2001), aimed at detect- 
ing "extended" infall in starless cores through the 
comparison of the line profiles of optically thick 
(CS) and thin (N 2 H + ) tracers. We considered all 
cores in the Lee, Myers, & Tafalla (2001) sam- 
ple showing strong N 2 H + (l-0) lines and compact 
morphology. Among these, 13 show infall asym- 
metry (i.e., CS spectra peaking at lower velocities 
than N 2 H + ), 7 show no asymmetry and 4 show 
outflow asymmetry. 

Of the other 7 cores (L1495, L1400K, TMC1C, 
L1689A, B68, L328, GF5) in our sample some have 
been mapped in lines of NH3 (Benson & Myers 
1989) and N 2 H+ (Caselli, Benson, Myers & Tafalla 
2002) and all have been mapped in dust contin- 
uum emission and show compact morphology and 
relatively high central number densities (n(H 2 ) > 
10 5 cm~ 3 ; Bacmann et al. 2000; Shirley, Evans, 
Rawlings, & Gregersen 2000; Tafalla et al. 2002). 

Finally, the cores in our sample are nearby 
enough (within 250 pc), to spatially resolve the 
dense core nucleus. The total sample is listed in 
Table 1; columns (2), (3) and (4) show the J2000.0 
coordinates and the distance in parsec, taken from 
the reference in column (5). In column (6) we 
specify if the core is isolated or embedded in a 
particular star forming region. Then in column 
(7) we give the total number of positions observed 
in both N 2 H+ and N 2 D+. Finally, in column (8) 
and (9) we report on the availability of a 1.2-mm 
map and of a C 18 map with the reference to it. 

2.2. N 2 H+ and N 2 D+ 

The N 2 H+ and N 2 D+ observations were carried 
out between April 2002 and January 2003 using 
the IRAM 30-m telescope. We observed in fre- 
quency switch mode and used the autocorrelator 
as the back end. Details on the telescope setting 
used are given in Table 2. 

The frequencies of N 2 H+(1 0), N 2 D+(2-l) 
and N 2 D + (3-2) were updated following the re- 
cent results of Dore et al. (2004). A new value 
for the N 2 H+(3-2) frequency was calculated, 
first evaluating the N 2 H + rotational constant (B 
= 46586.87546 ± 0.00025 MHz) from the new 
N 2 H + (l-0) frequency determination (Dore et al. 
2004), and then using the centrifugal distortion 
constant from Verhoeve et al. (1990). In Table 2, 
we report the frequency of the FiF = 45— »34 



hyperfine component which has a relative intensity 
of 17.46%. The pointing and focus were checked 
approximately every two hours observing planets 
or quasars in continuum. The pointing correc- 
tions never exceeded 4". The intensity scale was 
converted to main beam temperature according to 
the efficiencies reported on the IRAM 30-m Web 
site. 1 

We first observed the total sample at the peak 
position given either from previous N 2 H + ob- 
servations (Lee, Myers, & Tafalla 2001) or from 
the dust maps (Bacmann et al. 2000). Then, we 
made a map of the cores that showed the strongest 
N 2 D+(2-l) lines: L1521F, L1544, TMC2, L183, 
OphD, L492, L429 and L694-2. TMC2 was 
mapped only in N 2 D + and we thus evaluate 
the ratio of N 2 H + and N 2 D + towards the only 
available position observed in N 2 H+ (25" off the 
N 2 D+ strongest position). The N 2 H+ and N 2 D+ 
maps towards L1495, L1517B, L1498, L1400K and 
L1512 were taken as part of a different project on 
starless cores aimed at determining the ionization 
degree (P. Caselli et al., in preparation, T. Gatti 
et al., in preparation); here we just use the N 2 H+ 
and N 2 D + values at the dust peak (determined by 
Tafalla et al. 2002). The L1544 and L1521F data 
used in the present paper are taken from Caselli et 
al. (2002a) and Crapsi et al. (2004). respectively. 

2.3. Dust continuum and C ls O 

Millimeter continuum observations, taken with 
MAMBO (the IRAM 30-m bolometer), were al- 
ready available in the literature for 13 cores in our 
sample (L1544, OphD: Ward-Thompson, Motte 
& Andre 1999; L1689B, L328, L429, GF5: Bac- 
mann et al. 2000; L1498, L1495, L1517B, L1400K: 
Tafalla et al. 2002; L183: Pagani et al. 2003; 
L1521F: Crapsi et al. 2004; L694-2: M. Tafalla 
et al. 2005, in preparation). Moreover, a 1.2- 
mm map of B68 taken with SIMBA (the facility 
bolometer at SEST) was published in Bianchi et 
al. (2003). We obtained a 1.2-mm map of two more 
cores: TMC2, which presents high N 2 D + cmission 
and L492 which has very strong N 2 H + (l-0) lines. 
Observations were carried out by the IRAM staff 
between January 2003 and February 2003 using 
the new 117-channel bolometer. The IRAM 30- 
m half-power beamwidth (HPBW) at 1.2 mm is 

1 Available at http://www.iram.fr/IRAMES/index.htm 
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Table 2 

Telescope settings and parameters. 



line 


frequency 


HPBW 


-F ''throw 




res 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


N 2 H+(l-0) 
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0.042 
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16 
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0.038 
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10 
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0.050 


C 18 O(l-0) 


109.782173 


22 
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0.026 


C ls O(2-l) 


219.560357 


11 




450 


0.033 



Note. — Col. (2) line rest frequency (GHz); Col. (3) Half Power 
Beam Width (") Col. (4) Freq. Throw (kHz); Col. (5) System Temper- 
ature (K); Col. (6) Channel Spacing (km s _1 ) 



approximately 11". The calibration error is esti- 
mated to be ^10%. 

We mapped TMC2 in a 7'x7' region scanning in 
the azimuthal direction at 5"s _1 with a spacing 
between the scans of 9", reaching a sensitivity of 
6 mjy/beam at the map center. The zenith opti- 
cal depth, measured before and after the map, was 
0.2. L492 was mapped in a 3'x3' region scanning 
in the azimuthal direction at 5"s _1 with a spacing 
between the scans of 8"; the RMS at map cen- 
ter is 5 mjy beam -1 . The zenith optical depth, 
measured before taking the map, was 0.2. The 
data reduction (baseline fitting, sky-noise subtrac- 
tion and regridding) was performed using the NIC 
program of the GAG software developed at IRAM 
and the Observatoire de Grenoble (Broguie're et 
al. 2003). 

In the same fashion, we observed the C 18 O(l-0) 
and (2-1) for those cores where data were not al- 
ready available in literature, i.e., L183, L694-2, 
L492 and TMC2. We mapped in the on-the-fly 
mode in a region 150" x 150" wide (200" x 300" for 
L183 and 200" x 200" for TMC2) scanning in the 
azimuthal direction. Both the separation between 
two scans and the separation between two dumps 
were 5". C 18 data were already available for 
L1544 (Caselli et al. 1999), OphD, L328, L429 
(Bacmann et al. 2002), L1495, L1498, L1517B 
(Tafalla et al. 2002), B68 (Bergin, Alves, Huard, 
& Lada 2002) and L1521F (Crapsi et al. 2004). 



3. Results 

3.1. Spectra at dust peak 

In Figure 1 and 2, we show the observed spec- 
tra of N 2 H+(l-0) and N 2 D+(2-l) towards the 
core dust peaks. The cores where no line was 
detected were omitted. The data were reduced 
using CLASS, the line analysis program of the 
GAG software (Buisson et al. 2002). Notable 
in particular is the double-peaked profile with an 
infall asymmetry observed in L1544, as well as the 
red shoulder or wing seen in several other sources 
(OphD, L694-2, L492, L1521F, TMC1-C). In gen- 
eral, these more "complex" profiles are more ap- 
parent in the sources where the N 2 D + (2-l) line is 
strong (see below), although there are exceptions 
(e.g. TMC1 which has a relatively broad N 2 H + 
line with "outflow asymmetry" but weak N 2 D + (2- 
l)) 2 . The integrated intensities of these spectra 
plus the N 2 H+(3-2) and N 2 D+(3-2) are reported 
in Table 3. We also give there the offset of the 
N 2 D + (2-l) peak from the reference position given 
in Table 1. In our sample of 31 cores, we detected 
N 2 H+(l-0) emission with S/N>10 in 29 cores. 
LI 689 A and GF5, which have a weak detection, 
will not be considered further in this work. We de- 
tected N 2 D + (2-l) lines with signal-to- noise-ratio 
(S/N) higher than 10 in 10 cores (L1521F, TMC2, 
L1517B, L1544, L183, OphD, L492, L328, L429, 
L694); another 13 were detected with S/N higher 

2 The TMC1 position observed is close to the "CS peak" , sec 
Pratap et al. (1997). 
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than 4. In our calculations, as well as the error due 
to the RMS noise (shown in Table 3), we adopt a 
calibration uncertainty of 10% for N2H + (l-0) and 
15% for N 2 D+(2-l). 

From hyperfine structure fitting of the spectra 
in Figure 1 and 2, we derived the line parameters 
shown in Table 4. We report the observed LSR ve- 
locities and line widths in the two species and find 
in first approximation good agreement, suggesting 
that we are sampling the same regions along the 
line of sight. 

Total optical depth estimates are derived from 
a simultaneous fit to the hyperfine satellites as- 
suming that hyperfine levels are populated pro- 
portional to their statistical weights. These re- 
sults show that the total optical depth (i.e., the 
sum of the optical depths of the seven satellites) 
in N 2 H+(l-0) can reach values as high as 20; con- 
sidering the hyperfine structure of N 2 H + (1 0), a 
total optical depth of 20 means that the main com- 
ponent (F\ F = 2 3 — > 12), which has a fraction 
of 7/27 (= 25.9%) of the total line strength, is 
often optically thick while the weakest satellite 
(Fi F = 1 — > 11), with an intrinsic line strength 
of 1/27 (or 3.7%) of the total transition strength, 
is in almost all cases optically thin (its maximum 
opacity in our sample is 0.7, towards L183, L1521F 
and B68). We note, however, that the combina- 
tion of excitation anomalies (non-LTE effects in 
the hyperfine level populations, see Caselli, Myers 
& Thaddeus 1995) and highly non-Gaussian line 
profiles (e.g., L1544, B68, L492, L1521F, L694-2) 
increases the uncertainty in the opacity estimate 
from the simultaneous fitting of the hyperfine com- 
ponents. In fact, the fitting model adopted as- 
sumes constant abundance and excitation temper- 
ature along the line of sight and does not consider 
variations in velocity and line width. 

To crudely estimate the variations on the more 
complex spectra, we derived the opacity of the 
N 2 H+(l-0) spectra in B68 and L1544 by the fit 
of only a couple of hypcrfines and repeating that 
for all of the different couples. In this way, we 
found a scatter of a factor of 2 in the resulting 
opacities around the value derived by the simulta- 
neous fit of all seven components. The results of 
the hyperfine structure fitting of the high-S/N (> 
5) 1mm spectra are reported in Table 5. We also 
show in Figure 3 the C 18 O(l-0) and (2-1) spectra 
towards the peak of L1521F, L694-2, L183, TMC2 
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Fig. 1. — N 2 H + (l-0) spectra at peak sorted by de- 
creasing N 2 H + column density. The two sources 
with no detection (L1689A and GF5) were omit- 
ted. 
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Table 3 

Observed line integrated intensities at N 2 D+peak. 



core 


RA off 


DEC off 


I(N 2 H+(l-0)) 


I(N 2 H+(3-2)) 


I(N 2 D+(2-l)) 


I(N 2 D+(3-2)) 




// 


// 


K km s _1 


K km s — 1 


K km s~ 1 


K km s -1 


L1498 








2.27±0.01 


0.11±0.02 


0.12±0.02 


<0.01 


L1495 


40 


40 


3.18±0.03 


0.24±0.04 


0.23±0.05 


0.08±0.02 


L1495B 








1.35±0.05 


<0.09 


U.14±U.U4 


0.30±0.09 


L1495A-N 








3.4b±0.05 


<0.09 


0.18±0.02 


<0.09 


L1495A-S 








1.08±0.0o 


<0.12 


0.15±0.04 


<0.11 


L1521F 


-10 





5.8o±0.02 


0.74±0.04 


0.97±0.02 


0.29±0.01 


L1400K 


-20 


-40 


1.39±0.02 


0.10±0.04 


<0.04 


<0.02 


L 1400 A 








1.48±0.04 


0.25±0.12 


0.08±0.02 


0.13±0.07 


TMC2 





40 


3.48±0.03 


0.32±0.11 


0.66±0.04 


0.17±0.05 


TMC1 








3.51±0.06 


0.35±0.11 


0.20±0.03 


<0.10 


TMC1C 





40 


1.87±0.03 


<0.19 


0.24±0.04 


0.08±0.06 


L1507A 








1.78±0.05 


0.21±0.12 


0.15±0.03 


0.15±0.07 


CB23 








2.22±0.06 


0.20±0.14 


0.21±0.04 


0.37±0.10 


L1517B 


-15 


-15 


2.57±0.01 


0.11±0.02 


0.26±0.02 


<0.01 


L1512 








2.11±0.02 


<0.05 


0.17±0.03 


<0.02 


L1544 


20 


-20 


5.46±0.04 


0.56±0.09 


2.23±0.05 


0.55±0.03 


L134A 








1.13±0.03 


<0.25 


<0.04 


<0.07 


L183 


30 





4.20±0.01 


0.56±0.04 


1.63±0.02 


0.48±0.01 


OphD 


-20 


-50 


3.71±0.02 


0.39±0.05 


1.75±0.02 


0.59±0.02 


L1689A 








0.15±0.07 


<0.13 


<0.05 


<0.06 


L1689B 








1.27±0.05 


<0.13 


0.09±0.04 


0.12±0.06 


L158 








0.61±0.05 


<0.11 


<0.04 


<0.05 


L234E-S 








2.06±0.07 


0.40±0.26 


0.26±0.04 


<0.32 


B68 








2.37±0.01 


0.17±0.02 


0.13±0.03 


<0.03 


L492 


20 


20 


3.95±0.02 


<0.10 


0.30±0.02 


<0.04 


L328 








4.46±0.03 


0.42±0.09 


0.38±0.03 


0.19±0.04 


L429 


-20 


20 


6.11±0.02 


0.57±0.05 


2.03±0.02 


0.62±0.02 


GF5 








0.39±0.05 


<0.07 


<0.04 


<0.10 


L694-2 








5.13±0.01 


0.36±0.02 


1.35±0.02 


0.37±0.01 


L1197 








2.15±0.04 


<0.16 


0.18±0.06 


<0.05 


CB246 


-60 


-20 


2.91±0.03 


<0.07 


0.11±0.04 


0.22±0.06 



Note. — Offsets are referred to the positions in Table 1. Temperatures are in main beam temperature scale. 
Calibration errors were not included in the error estimate. Upper limits assume the line-width of N2H+(l-0) 
where available or 0.17 km s _1 the smallest line-width observed in our sample. 
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Table 4 

Results of the hyperfine structure fitting on the N 2 H + (l-0) and N 2 D+(2-l) spectra at 

THE PEAK POSITION OF EACH CORE. 







km s _1 




N 2 H+(l-0) 
AV t 

km s -1 




i ■ 

1 EX 

K 






km s _1 




N 2 D+(2-l) 
AV 
km s _1 




^ EX 

K 




L1498 


7 


822±0 


001 





185±0 


001 


11.5±0.3 


4 


5±0 


1 


7 


778±0 


010 





176±0 


029 




0.1 


4.5 




L1495 


6 


807±0 


001 





237±0 


003 


9.6±0.6 


4 


8±0 


2 


6 


755±0 


025 





290±0 


061 




0.1 


4.8 




L1495B 


6 


633±0 


008 





390±0 


018 


0.1 




4.8 




6 


683±0 


065 





500±0 


216 




0.1 


4.8 




L1495A-N 


7 


296±0 


003 





327±0 


007 


7.6±0.9 


4 


7±0 


3 


7 


237±0 


013 





244±0 


035 




0.1 


4.7 




L1495A-S 


7 


294±0 


007 





214±0 


015 


21.3±5.3 


4 


4±0 


4 


7 


264±0 


024 





284±0 


059 




0.1 


4.4 




L1521F 


6 


472±0 


001 





299±0 


001 


17.9±0.3 


4 


9±0 


1 


6 


505±0 


004 





268±0 


010 


2 


2±0.4 


4.6±0 


3 


L1400K 


3 


196±0 


002 





231±0 


005 


6.5±0.8 


4 


4±0 


5 






















L1400A 


3 


355±0 


002 





191±0 


006 


7.7±1.2 


4 


4±0 


6 


3 


250±0 


036 





175±0 


061 




0.1 


4.4 




TMC2 


6 


193±0 


001 





210±0 


002 


8.7±0.4 


5 


1±0 


2 


6 


157±0 


006 





167±0 


015 


4 


5±1.1 


5.0±0 


3 


TMC1 


5 


856±0 


003 





269±0 


007 


11.1±1.2 


4 


6±0 


3 


5 


870±0 


040 





409±0 


113 




0.1 


4.6 




TMC1C 


5 


196±0 


003 





212±0 


006 


18.4±2.2 


4 


4±0 


5 


5 


193±0 


009 





132±0 


017 




0.1 


4.4 




L1507A 


6 


163±0 


004 





220±0 


009 


12.2±2.1 


4 


4±0 


6 


6 


159±0 


030 





254±0 


071 




0.1 


4.4 




CB23 


6 


015±0 


002 





166±0 


006 


14.6±2.2 


4 


4±0 


3 


5 


985±0 


027 





233±0 


075 




0.1 


4.4 




L1517B 


5 


835±0 


001 





215±0 


001 


8.8±0.2 


4 


7±0 


1 


5 


796±0 


008 





218±0 


017 




0.1 


4.7 




L1512 


7 


121±0 


001 





174±0 


001 


7.3±0.4 


5 


0±0 


2 


7 


091±0 


018 





274±0 


076 




0.1 


4.8 




L1544 


7 


143±0 


002 





315±0 


004 


12.6±0.7 


5 


0±0 


2 


7 


181±0 


004 





289±0 


008 


5 


2±0.4 


4.8±0 


2 


L134A 


2 


665±0 


002 





215±0 


007 


4.3±1.0 


4 


4±0 


6 






















L183 


2 


413±0 


001 





211±0 


001 


20.3±0.2 


4 


8±0 


1 


2 


459±0 


001 





214±0 


002 


4 


7±0.2 


4.9±0 


1 


OphD 


3 


478±0 


001 





223±0 


001 


4.4±0.1 


7 


1±0 


2 


3 


515±0 


001 





204±0 


003 


4 


3±0.2 


5.3±0 


2 


L158 


3 


942±0 


007 





229±0 


014 


0.1 




4.8 
























L1689B 


3 


481±0 


005 





262±0 


011 


0.1 




4.8 




3 


642±0 


058 





373±0 


148 




0.1 


4.8 




L234E-S 


3 


164±0 


003 





198±0 


008 


6.9±1.5 


4 


7±0 


5 


3 


286±0 


015 





182±0 


055 




0.1 


4.7 




B68 


3 


364±0 


001 





176±0 


001 


26.2±1.0 


4 


5±0 


5 


3 


372±0 


016 





124±0 


036 




0.1 


4.5 




L492 


7 


701±0 


001 





263±0 


001 


11.8±0.2 


4 


8±0 


1 


7 


726±0 


008 





222±0 


022 


2 


6±0.8 


4.4±0 


3 


L328 


6 


707±0 


002 





438±0 


005 


8.1±0.5 


4 


6±0 


2 


6 


797±0 


012 





419±0 


028 




0.1 


4.6 




L429 


6 


719±0 


001 





394±0 


001 


12.3±0.1 


4 


9±0 


1 


6 


705±0 


002 





360±0 


005 


4 


3±0.2 


4.4±0 


1 


L694-2 


9 


574±0 


001 





266±0 


001 


13.9±0.1 


5 


2±0 


1 


9 


567±0 


002 





244±0 


004 


5 


5±0.2 


4.1±0 


1 


L1197 


-3 


.147±0 


.002 





249±0 


006 


6.5±0.9 


4 


5±0 


3 


-3 


.020±0 


.016 





180±0 


093 




0.1 


4.5 




CB246 


-0 


.830±0 


.001 





245±0 


003 


10.3±0.7 


4 


5±0 


2 


-0 


.730±0 


.025 





126±0 


039 




0.1 


4.5 





Note. — The hyperfine fitting model assumes constant Tex and homogenous abundance and do not consider velocity 
gradients effects along the line of sight. L1689A and GF5 spectra were not fitted because of poor S/N. 
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Table 5 

Results of the hyperfine structure fitting on the N 2 H+(3-2) and N 2 D+(3-2) spectra at 

THE PEAK POSITION OF EACH CORE. ONLY THE FITS TO HIGH S/N (> 5) SPECTRA ARE REPORTED. 







N 2 H+(3-2) 






N 2 D+(3-2) 




core 


VlS-R 




T 




AV 


T 




km s _1 


km s _1 




km s — 1 


km s — 1 




L1498 


7.786±0.015 


0.169±0.040 


<0.1 








L1495 


6.811±0.030 


0.288±0.052 


<0.1 








L1521F 


6.391±0.010 


0.172±0.021 


6.2±1.5 


6.507±0.006 


0.222±0.025 


1.2±0.9 


L1517B 


5.816±0.030 


0.219±0.130 


<0.2 








L1544 


7.138±0.018 


0.193±0.080 


<0.1 


7.208±0.012 


0.342±0.026 


<0.1 


L183 


2.412±0.012 


0.220±0.036 


3.5±1.6 


2.445±0.003 


0.205±0.012 


1.2±0.4 


OphD 


3.466±0.011 


0.131±0.019 


6.7±2.2 


3.503±0.004 


0.196±0.013 


2.6±0.6 


B68 


3.370±0.017 


0.231±0.022 


0.8±0.4 








L429 


6.673±0.029 


0.373±0.037 


2.7±0.4 


6.696±0.009 


0.262±0.020 


5.2±0.9 


L694-2 


9.532±0.014 


0.308±0.060 


1.7±1.4 


9.560±0.005 


0.242±0.016 


1.5±0.5 



Table 6 

Integrated intensities and results of the gaussian fitting of the C 18 O(l-0) and (2-1) 

SPECTRA TOWARDS THE PEAK POSITION FOR THE FOUR CORES NEWLY OBSERVED. 







C 18 O(l-0) 






C ls O(2-l) 




core 


I(C 18 O(l-0)) 




AV 


I(C ls O(2-l)) 


Vlsr 


AV 




K km s _1 


km s~ 1 


km s^ 1 


K km s _1 


km s~ 1 


km s~ 1 


L1521F 


1.80±0.06 


6.458±0.006 


0.438±0.014 


1.80±0.15 


6.427±0.015 


0.447±0.034 


TMC2 


1.05±0.04 


6.313±0.007 


0.401±0.014 


1.09±0.07 


6.229±0.010 


0.392±0.023 


L183 


1.89±0.04 


2.518±0.006 


0.545±0.016 


2.14±0.06 


2.556±0.008 


0.579±0.019 


L492 


0.76±0.04 


7.786±0.017 


0.651±0.042 


0.58±0.05 


7.786±0.022 


0.529±0.055 


L694-2 


1.21±0.04 


9.574±0.007 


0.433±0.015 


1.01±0.04 


9.574±0.006 


0.415±0.015 
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L492 n 

L1517B .v^Vjv^, 

CB23 

TMC1 



L1495A-N ^. 
L1507A 
L1 495A-S 
L1512 
B68 



Velocity (km/s) 



Fig. 2. — N 2 D + (2-l) spectra at peak sorted by de- 
creasing N2D" 1 " column density. Eight sources with 
marginal or no detection were omitted (L1400K, 
L1400A, L134A, L1689A, L1689B, L158, GF5, 
CB246). 




-2 2 

Velocity (km/s) 

Fig. 3.— C 18 O(l-0) and (2-1) spectra towards 
L1521F, L694-2, L183, TMC2 and L492. The 
shown spectra are the arithmetic means of all the 
observations in the 10" strip centered around the 
peak position and going parallel to the Right As- 
cension direction. 
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and L492. Integrated intensities and results of the 
Gaussian fits for these observations are reported 
in Table 6. 

3.2. Maps of individual sources 

A detailed description of the maps obtained in 
this work can be found in the Appendix A, where 
we also analyse core shapes and kinematics (Ap- 
pendix B) of individual objects. Here we briefly 
summarize the main results. Our new maps con- 
firm the well known finding that N2H" 1 " and N2D + 
well trace the millimeter dust continuum emission, 
unlike CO (e.g. Caselli et al. 1999; Bergin et al. 
2001; Tafalla et al. 2002). We also find that the 
majority of the cores are "cometary" in shape, sug- 
gesting that the cores have been either compressed 
by external agents or have moved relatively to the 
surrounding environment (but see Walsh & Myers 
2004). In general, they do not show an axisymmet- 
ric morphology. Finally, we studied the kinematics 
in individual cores, evaluating the "total" and "lo- 
cal" velocity gradients and finding that only in one 
case (L694-2) the N2H + (l-0) line width increases 
toward the core center, indicating the presence of 
enhanced dynamical activity (likely due to infall) , 
as already found in L1544 (Caselli et al. 2002a) 
and L1521F (Crapsi et al. 2004). 

In the rest of the paper we concentrate on the 
physical and chemical parameters at the relative 
peak of the whole source sample. 

4. Analysis 

In this section we derive the physical and chemi- 
cal properties that can best discriminate the evolu- 
tionary status of a starless core, either in a chem- 
ical sense, as a more pronounced deuterium en- 
richment and a higher CO depletion factor, or for 
dynamical reasons, like the presence of a denser 
and more centrally concentrated H2 distribution 
or broader N 2 H + line widths. We expect to char- 
acterize through these parameters the pre-stellar 
cores, i.e., the subset of starless cores that present 
enhanced chemical evolution and/or dynamical ac- 
tivity. 

4.1. Column densities of N 2 H+ and N 2 D+ 
and the deuterium fractionation 

The deuterium fractionation can be evaluated 
from observations of a hydrogen bearing molecule 




5 10 15 20 

N(N 2 H+) [10 12 cm- 2 ] 



Fig. 4. — N 2 H + column density vs. N 2 D+ col- 
umn density at N2D + peak position. Deuterium 
fractionation can be read from the radial dashed 
lines (corresponding to 7V(N2D+)/7V(N2H+) of 
0.1, 0.25 and 0.5). Dotted lines are the arithmetic 
means of the distributions. Cores belonging to the 
Taurus Molecular Cloud are marked by a square. 



and one of its deuterated counterparts by calcu- 
lating the ratio of their column densities. Re- 
cent work has made use of HCO + , H2CO, HNC, 
N 2 H+ and NH 3 , as well as their single or multiple 
deuterated isotopologues, to evaluate this param- 
eter in starless cores (Caselli, Walmsley, Terzieva 
& Hcrbst 1998; Roueffet al. 2000; Tine et al. 2000; 
Bacmann et al. 2003; Caselli et al. 2002b; Hirota, 
Ikeda, & Yamamoto 2003; Crapsi et al. 2004). 
One should note, however, given the depletion of 
HCO + and H 2 CO in the core nucleus (Caselli et 
al. 2002b; Carey et al. 1998), that the deuterium 
fractionation estimated from this species is rela- 
tive to an outer shell and not to the high-density 
nucleus. 

In this work we derived N 2 H+ and N 2 D + col- 
umn densities assuming a constant excitation tem- 
perature (the CTEX method: see the appendix 
of Caselli et al. 2002b). The excitation temper- 
ature of N 2 H + (l-0) was derived from hyperfinc 
fitting of the spectra with well constrained opaci- 
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ties (r/(j T > 3, where the total optical depth r and 
its associated error a T have been derived assuming 
a constant excitation temperature for the hyper- 
fine components; see also Table 4); otherwise the 
spectrum was assumed optically thin and the exci- 
tation temperature was set to 4.8 K (the average 
excitation temperature in our sample). The ex- 
citation temperature of N 2 D + (2-l) was set equal 
to that of N 2 H + (l-0) in those cases in which the 
opacity was not well constrained. The effect of un- 
derestimating the N 2 D + (2-l) excitation tempera- 
ture by, e.g., 1 K would be a decrease of 7V(N 2 D+) 
by 30%. 

As noted above, in the case of N 2 H + , the optically 
thin assumption is often incorrect and hence we 
have used the intensity of the (F\ F = 1 — > 11) 
(weak) component normalized appropriately in or- 
der to estimate the J=0 column density. 

The final result can be seen in Figure 4, 
where the N 2 H + column density is plotted ver- 
sus the N 2 D+ column density. The estimated 
V(N 2 D+)/7V(N 2 H+) ratio varies between 0.44 in 
Oph D and an upper limit of 0.02 in L1400K. It 
is noticeable that with the exception of OphD, all 
cores with 7V(N 2 D+)/7V(N 2 H+) ratio above 0.10 
have N 2 H + column densities above 10 13 cm~ 2 , 
whereas cores with N (N 2 D+)/7V(N 2 H+) lower 
than 0.10 have V(N 2 H+) lower than 10 13 cm" 2 
(with a few exceptions). Thus, the impression 
from Figure 4 is that there is a trend in the sense 
that the V(N 2 D+)/7V(N 2 H+) ratio correlates pos- 
itively with N 2 H + column density. Table 7 sum- 
marize the above numbers. 

We also determined the N 2 H + and N 2 D + col- 
umn densities in the Large Velocity Gradient 
(LVG) approximation in those cores with high- 
S/N N 2 H+ and N 2 D+(3-2) spectra (see Sec- 
tion C). Results are shown in Table 12. The 
LVG column density determinations are found to 
agree with the CTEX values with a dispersion of 
30%. 

4.2. The integrated CO depletion factor 

The CO integrated depletion factor is the ra- 
tio of the CO canonical abundance ([CO]/[H 2 ] = 
9.5 x 10~ 5 , Frerking, Langcr, & Wilson 1982) and 
the CO observed abundance. The CO observed 
abundance integrated along the line of sight can 
be calculated from the ratio of the CO column 
density and the H 2 column density. C 18 is mea- 



sured to be optically thin in some of the cores 
(L1521F, L183, L694-2, TMC2, L492) where C 17 
observations are available. Optical depth was de- 
rived from the integrated intensity ratio of the 
two isotopologues assuming equal excitation tem- 
perature and the local interstellar medium rela- 
tive abundance value (3.2, Wilson & Rood 1994). 
The maximum optical depth was measured to be 
0.8 towards L1521F. Thus, for practical purposes, 
we assumed C 18 optically thin in all the cores. 
To derive V(CO), we used the CTEX approxi- 
mation on the C 18 lines. As already shown in 
Caselli et al. (2002b) and Crapsi et al. (2004), 
this approach reduces the evaluation of the inte- 
grated CO depletion factor to the simple formula 
f D = 0.085 • T x . 2mm/^Vc 18 Oi where -7"i. 2 mm is the 
1.2-mm observed brightness in mJy beam -1 and 
Wqwq is the C 18 integrated intensity in K km 
s- 1 . 

The sources of uncertainties in this technique 
are several, and it is important to keep them in 
mind in the analysis of the results. First, the 
"canonical CO abundance" measured by different 
authors or in different objects varies by a factor of 
2 (Lacy et al. 1994; Alves, Lada, & Lada 1999; 
Kramer et al. 1999). Second, one has to bear in 
mind that the integrated CO depletion factor is an 
average value along the line of sight, thus the local 
depletion factor in the nucleus of the core is much 
greater (see, e.g., Caselli et al. 2002b; Crapsi et 
al. 2004). 

The integrated depletion factors of L1689B, 
L328, OphD and L429 were already published in 
Bacmann et al. (2002), that of L1544 in Caselli 
et al. (2002b) and that of L1521F in Crapsi et 
al. (2004). We evaluate the f D values for L1495, 
L1498 and L1517B using the data published in 
Tafalla et al. (2002), for B68 using the C 18 O(l-0) 
intensities published in Bergin, Alves, Huard, & 
Lada (2002) and the 1.2-mm continuum data pub- 
lished by Bianchi et al. (2003), for L183 using the 
1.2-mm map of Pagani et al. (2003) and our C 18 
observations, for L694-2 using the 1.2-mm map of 
M. Tafalla et al. (2005, in preparation) and our 
C 18 observations and for TMC2 and L492 using 
the C 18 and the 1.2-mm maps from this work. 
Results are shown in column (4) of Table 8. 
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Table 7 

n 2 h+ and n 2 d+ column densities and relative ratio towards the n 2 d + (2-l) integrated 

intensity peak. 



core 


iV(N 2 H+) 
10 12 cm- 2 


Ar(N 2 D+) 
10 11 cm- 2 


A r (N 2 D+)/jV(N 2 H+) 


L1498 


7.1±0.7 


2.7±0.6 


0.04±0.01 


L1495 


8.8±1.0 


4.2±1.1 


0.05±0.01 


L1495B 


2.8±0.8 


2.7±0.8 


0.10±0.04 


L1495A-N 


9.6±1.2 


3.6±0.7 


0.04±0.01 


L1495A-S 


4.6±1.0 


3.4±1.0 


0.08±0.03 


L1521F 


19.3±1.9 


20.0±3.0 


0.10±0.02 


L1400K 


3.6±0.5 


<0.8 


<0.02 


L 1400 A 


3.6±0.7 


1.8±0.6 


0.05±0.02 


TMC2 


9.9±1.1 


11.3±1.8 


0.11±0.02 


TMC1 


9.6±1.3 


4.1±0.9 


0.04±0.01 


TMC1C 


8.0±0.9 


5.5±1.3 


0.07±0.02 


L1507A 


6.9±0.9 


3.5±0.9 


0.05±0.01 


CB23 


7.0±1.2 


4.8±1.1 


0.07±0.02 


L1517B 


7.9±0.8 


5.1±0.8 


0.06±0.01 


L1512 


5.8±0.6 


3.1±0.6 


0.05±0.01 


L1544 a 


18.3±1.9 


42.5±6.4 


0.23±0.04 


L134A 


2.5±0.5 


<0.9 


<0.04 


L183 


13.7±1.4 


29.9±4.5 


0.22±0.04 


OphD 


6.3±0.7 


27.3±4.1 


0.44±0.08 


L158 


1.0±0.1 


<0.7 


<0.08 


L1689B 


2.0±0.2 


1.7±0.8 


0.09±0.04 


L234E-S 


5.9±1.1 


5.0±1.1 


0.08±0.02 


B68 


9.7±1.0 


2.7±0.7 


0.03±0.01 


L492 


14.4±1.5 


6.9±1.1 


0.05±0.01 


L328 


11.6±1.3 


7.9±1.3 


0.07±0.01 


L429 


16.7±1.7 


45.9±6.9 


0.28±0.05 


L694-2 


14.6±1.5 


37.6±5.7 


0.26±0.05 


L1197 


5.7±0.9 


4.0±1.5 


0.07±0.03 


CB246 


8.8±1.0 


2.4±0.9 


0.03±0.01 



a The difference in the L1544 N (N 2 D+)/iV (N 2 H+) given by Caselli et al. 
(2002b) (0.26) and here (0.23) is due to the different approximation used. They 
evaluated the column density from all the seven hyperfine components correcting 
for the high optical depth (see their eq. Al) while here we used the optically 
thin approximation applied on the weakest component only. 
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Table 8 

h 2 densities, integrated co depletion factor, flattened radius, aspect ratios and 

N 2 H+(l-0) SKEWNESS. 



core 


JV(H 2 ) 
10 22 cm- 2 


n(H 2 ) 
10 5 cm" 3 


f D (CO) a 


T-70 

10 3 AU 


aspect ratio 
70% contour 


aspect ratio 
50% contour 


N 2 H+(l-0) 
skewness b 


L1498 


3.2±1.0 


1.0±0.7 


7.5i2.5 


6.6i0.7 


1.6i0.3 


1.5i0.3 


-0.01i0.05 


L1495 


3.1±1.0 


1.1±0.7 


7i2.4 


4.7i0.6 


2.0i0.4 


1.9i0.4 


0.10i0.05 


L1521F 


13.5±2.2 


11. Oil. 8 


15i3.6 


3.4i0.5 


1.2i0.1 


1.6i0.2 


0.19i0.06 


TMC2 


6.0±1.2 


3.0i0.8 


13i3.3 


5.2i0.6 


1.2i0.4 


l.li0.3 


-0.01i0.05 


L1517B 


3.7±1.0 


2.2i0.8 


9.5i2.8 


4.0i0.5 


1.4i0.1 


1.2i0.1 


-0.08i0.05 


L1544 


9.4±1.6 


14.0±2.2 


14i3.4 


3.2i0.4 


1.6i0.2 


1.9i0.2 


0.40i0.09 


L183 


10.0±1.7 


10.0±1.7 


12i3.1 


4.8i0.6 


1.8i0.2 


2.5i0.2 


0.09i0.05 


OphD 


8.2±1.5 


8.5±1.5 


14i3.4 


6.1i0.7 


2.1i0.4 


1.9i0.4 


0.54i0.12 


B68 


1.4±0.3 


0.8i0.7 


3.4i2.1 


3.7i0.9 


1.3i0.3 


1.3i0.3 


-0.09i0.05 


L492 


4.4±1.1 


2.1±0.8 


8i2.6 


7.8i0.8 


1.5i0.2 


2.0i0.2 


0.25i0.07 


L328 


5.7±1.2 


1.8±0.7 


8.5i2.6 


3.7i0.5 


l.li0.2 


1.2i0.2 


-0.02i0.08 


L429 


8.8±1.6 


6.0±1.1 


15.5i3.7 


3.6i0.5 


1.4i0.2 


2.1i0.3 


-0.20i0.10 


L694-2 


7.8±1.4 


9. Oil. 5 


lli3.0 


5.5i0.6 


1.3i0.1 


1.4i0.1 


0.22i0.07 



Note. — For reference to the literature data see text. 
a The CO depletion factor in L1689B is 4.5i2.2. 

b N 2 H+(l-0) skewness was measured also in: L1495A-N (-0.04i0.05); L1400K (-0.42i0.10); TMC1 (- 
0.13i0.06); TMC1C (0.43i0.10); L1512 (0.02i0.05); L234E-S (-0.25i0.07); L1197 (0.08i0.11); CB246 
(0.02i0.05). 



4.3. H 2 volume density 

The H 2 column density can be derived from 
the millimeter continuum under the approxima- 
tion of optically thin emission, constant tempera- 
ture and constant emissivity of the dust at these 
wavelengths. The equation that relates the H 2 col- 
umn density to the 1.2-mm flux is: 



F>i,\± ) I-beam Kl.2mm 

where iV(H 2 ) is the H 2 column density, Si. 2 mm 
is the flux density integrated over the solid angle 
^beam) B V {T) is the Planck function at tempera- 
ture T, ni.2mm is the dust opacity per unit mass 
and m is the mean molecular mass. In our calcula- 
tions, we assume a dust temperature of 10 K and a 
dust opacity of 0.005 cm 2 g _1 . The assumed dust 
temperature is justified by the estimates obtained 
from NH 3 (Tafalla et al. 2002) and H 2 CO (Bac- 
mann et al. 2002) observations in several cores of 
our sample. We note that the dust temperature is 
predicted to drop in the core nucleus (Evans et al. 
2001; Zucconi et al. 2001); assuming a tempera- 
ture of 8 K would cause an increase of the column 



density of a factor of 1.5. We also remark that 
the dust opacity value, assumed following Andre, 
Ward-Thompson & Motte (1996), suffers from an 
uncertainty of a factor of 2 (Hcnning, Michel, & 
Stognienko 1995; Bianchi et al. 2003; Kramer et 
al. 2003). 

In this work we both collected 1.2-mm maps 
from the literature and observed new sources with 
the bolometers. In Table 8 we report the H 2 vol- 
ume density at the peak of L429 and L328 taken 
from Bacmann et al. (2000), of L1495, L1498, 
L1517B and L1544 from Tafalla et al. (2002), of 
OphD and L1689B from Ward-Thompson, Motte 
& Andre (1999), of L694-2 from Harvey et al. 
(2003) (using their Bonnor-Ebert fit and chang- 
ing Ki.2mm to 0.005 cm 2 g _1 for consistency), of 
L183 from Pagani et al. (2003) and of L1521F from 
Crapsi et al. (2004). All the central densities above 
were determined from 1.2-mm continuum data and 
using the same basic technique (i.e., a fit of the 
observed continuum emission starting from a vol- 
ume density model of the core). We remark that 
these estimates were obtained assuming spherical 
symmetry. A strong deviation of the density pro- 
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file along the line of sight with respect to what 
is seen in the plane of the sky will cause errors 
in the volume density determination. For some of 
these cores other different estimates for the peak 
density are reported in the literature; these were 
obtained starting from different data (molecular 
line or 850 /im continuum) or using different tech- 
niques. In the case of multiple estimates all done 
starting from 1.2- mm continuum data, we adopted 
the most recent value, e.g., L1544 has three deter- 
minations of the density that were done from a 
1.2-mm map: 1.5xl0 6 cm -3 by Ward-Thompson, 
Motte & Andre (1999), 4xl0 5 cm' 3 by Bacmann 
et al. (2000) and 1.4xl0 6 cm" 3 by Tafalla et al. 
(2002); in this work we adopted the most recent 
determination from Tafalla et al. (2002). 

The H 2 central volume density at the TMC2, 
L492 and B68 peaks was evaluated with the same 
technique as for L1521F (Crapsi et al. 2004). We 
assumed the analytical model for the volume den- 
sity given by Tafalla et al. (2002): n(H 2 ) = 
no/(l + (r/r ) a ), where the central density (n ), 
the "flattened radius" (r ) and the steepness of 
the profile at large radii (a) are free parameters. 
Then, we searched for the best combination of pa- 
rameters able to predict the observed continuum 
measurements. In this way we found at the peak 
of TMC2 an H 2 density equal to 3 x 10 5 cm" 3 , 
while 2.1 x 10 5 cm~ 3 is the measured density to- 
wards the L492 peak. The B68 central density, 
derived from lower S/N data (and with the worst 
angular resolution), was found to be 8 x 10 4 cm~ 3 
with large error bars. 

Volume densities were also derived using an 
LVG program (see Section C), whose results are 
summarized in Table 12. We found systemati- 
cally lower central densities using the molecular 
data (~ 3 times less) than using the dust (as al- 
ready found in L1544; see Caselli et al. 2002b), 
suggesting that N 2 H + and N 2 D+ may be probing 
a shell exterior to the high-density nucleus (see 
Figure 17 in Section C). This indication is con- 
sistent with the idea that N 2 could freeze-out for 
densities higher than 5 x 10 5 cm -3 as suggested by 
Bergin, Alves, Huard, & Lada (2002) and Bcllochc 
& Andre (2004). Keto, Rybicki, Bergin, & Plume 
(2004) found a similar result using N 2 H + and sug- 
gest that a way to reconcile gas density estimates 
from dust observations with those from molecu- 
lar line observations is to increase the dust mass 



opacity to k = 0.04 cm 2 g _1 , a value typical of 
fluffy aggregates with ice mantles. We stress that 
higher resolution observations are strongly needed 
to confirm this trend. 

4.4. Dust emission equivalent radius and 
aspect ratio 

In starless cores the H 2 density profile can be 
approximated by a region of roughly constant 
density followed by a power-law decrease (Ward- 
Thompson, Motte & Andre 1999; Bacmann et al. 
2000; Tafalla et al. 2002): this behaviour is re- 
flected in the 1.2-mm continuum emission. Since 
the "flattened nucleus" is expected to become 
smaller with the ongoing of infall, its size gives us 
a measure of the relative contraction stage reached 
by the cores. Unfortunately each author uses dif- 
ferent definitions for the size of the flattened re- 
gion; thus, we have measured, for the 13 cores with 
a 1.2-mm map, the area within the 70% contour 
of the dust peak. We then translated this into a 
radius (reported as rj n in figures and tables) by 
evaluating the square root of this area divided by 
7r. The resulting values are reported in column (5) 
of Table 8. This approach is useful for a relative 
comparison of the different cores. The choice of 
the 70% contour as threshold is justified by the 
fact that we wanted to limit the influence of the 
"cometary tail" (see Section B) on the r? deter- 
mination. 

We also evaluated the aspect ratio of the dust 
emission measuring the major and minor axis size 
of the 50% contour and 70% contour (results are 
shown in columns (6) and (7) of Table 8). We 
found that the axis ratio increases at lower con- 
tours in the more cometary cores, reaching values 
up to 2.5. 

4.5. The line asymmetry 

The analysis of the line profile gives us clues on 
the core kinematics. In fact, motions along the 
line of sight combined with optical depth effects, 
molecular abundance variations and temperature 
gradients alter the line profile, causing departures 
from a Gaussian line shape. In particular, looking 
at the density peak, we expect a "blue shoulder" in 
the case of outflow motions and a "red shoulder" 
in the case of infall motions (see Evans 2002, for 
a review). 
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Here, we determine the degree of line asymme- 
try by calculating the line skewness, or the third 
moment, of a distribution equivalent to the line 
profile. This is defined as X^ili (( x i ~ S)/ ^) 3 /N, 
where x and a x are the first and second moments 
of the same distribution and N is the total num- 
ber of points (see e.g. Press, Teukolsky, Vetter- 
ling & Flannery 1992). By definition, the third 
moment measures the deviation from a symmet- 
ric distribution, as the first moment furnishes its 
mean and the second yields its dispersion. It is 
a non-dimensional quantity and assumes positive 
values for distributions that have a "tail" in the 
values higher than the mean and negative values 
in the opposite case. We applied this analysis to 
the (F 1 F = 2 2 -> 1 1) component of N 2 H+ (the 
fourth starting from the lower velocities). The 
"isolated" (F 1 F = 01 -> 12, the one at low- 
est velocity) and the "weak" (F\ F = 1 — > 11, 
the one at highest velocity) components show ap- 
proximately the same behaviour as this one in the 
very high S/N spectra (RMS < 0.03 K with a 
0.063 km s _1 channel spacing); thus, in order to 
enlarge the sample to spectra with slightly lower 
S/N (RMS < 0.1 K), we used the data from the 
(iq F = 2 2 — > 1 1) component, which also pro- 
vides a sufficient isolation from the other com- 
ponents even in the spectra with broadest line 
widths. Moreover, given the role of the optical 
depth in the production of the asymmetric pro- 
files, the (F\ F = 2 2 — > 1 1) component is sup- 
posed to have a more pronounced skewness than 
the other two components. Results are shown in 
Table 8. Eight cores show red-skewed "contract- 
ing" profiles, six cores show blue-skewed "expand- 
ing" profiles, and seven cores have symmetric pro- 
files within the error bars. Note that the line peak 
in the "contracting" profiles falls at lower veloci- 
ties (or on the blue side) than x whereas asymmet- 
ric profiles with outflow character are the reverse. 

5. Discussion 

In the following we discuss the relations be- 
tween the quantities determined in the previous 
sections and compare them with simple theoreti- 
cal expectations. We divide this section into two 
main parts: in the first we concentrate on the con- 
nection between deuterium fractionation and the 
other chemical and dynamical parameters derived 
in Section 4, while in the second we cross-correlate 



only the parameters linked with the kinematical 
activity. 

5.1. Deuterium Fractionation correlations 

5.1.1. Deuterium fractionation and CO deple- 
tion 

An observational proof of the chemical model 
relating the deuterium fractionation with the CO 
depletion factor was given in Bacmann ct al. 
(2003) using the ratio of D 2 CO over H 2 CO in a 
sample of five cores. Although D 2 CO and H 2 CO 
have been shown to deplete in the inner nucleus 
of the pre-stellar cores (Carey et al. 1998; Maret 
ct al. 2004) and hence do not trace the very inner 
nucleus, these five cores (included in our sample) 
showed a good correlation between the CO de- 
pletion and the deuterium fractionation in H 2 CO 
with the exception of OphD. 

Another possible test is to check for a corre- 
lation between the CO depletion factor and the 
deuterium fractionation within a given core. This 
test was performed in Crapsi et al. (2004) showing 
again an affirmative answer. 

Here we searched for a correlation between the 
integrated CO depletion factor and the deuterium 
fractionation in a subsample of 14 cores using 
the high density gas tracers N 2 H + and N 2 D + 
to derive the deuterium fractionation. The re- 
sult is shown in Figure 5, where we denote by 
enclosed data points our "candidate pre-stellar 
cores" , defined for the purpose of this work as hav- 
ing iV(N 2 D+)/7V(N 2 H+) > 0.1. 

The correlation is not extremely tight (corre- 
lation coefficient 69%; 74% if only Taurus cores 
are considered), but indeed we can say that the 
cores that show higher CO depletions do also show 
higher deuterium fractionation. As in Bacmann et 
al. (2003), OphD seems to have "too much" deu- 
terium fractionation compared to the CO deple- 
tion. 

In Figure 5 we also show theoretical curves from 
simple chemical models based on that described 
by Caselli et al. (2002b). These models all assume 
that the gas and dust temperatures are constant 
across the core and equal to 10 K. The chemi- 
cal network and chemical parameters are also the 
same in all the models: apart from H 2 , the chemi- 
cal network contains the three neutral species CO, 
N 2 , and O, which can freeze out onto dust grains 
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Fig. 5. — Deuterium fractionation vs. integrated 
CO depletion factor at the N 2 D + peak position. 
Depletion factor values were taken from the liter- 
ature and from present data (see Section 4.2 for 
references). The solid line connects the predic- 
tions of chemical models of spherically symmetric 
dense cores with various degrees of central concen- 
trations (from central densities of 4.4 x 10 4 cm' 3 in 
model m to 4.4xl0 8 cm~ 3 in model n 5 ), roughly 
simulating (in spherical symmetry) the evolution- 
ary sequence of contracting disk-like clouds in the 
Ciolek & Basu (2000) model (see text). Dotted 
lines are the arithmetic means of the distributions. 
Cores belonging to the Taurus Molecular Cloud 
are marked by a square; pre-stellar cores candi- 
dates are circled. 



and return to the gas phase via thermal desorp- 
tion or cosmic- ray impulsive heating (following 
the formulation of Hasegawa & Herbst 1993). On 
the other hand, the abundances of the molecular 
ions (N 2 H+, HCO+, H 3 0+, H+ and their (mul- 
tiply) deuterated counterparts ) are calculated in 
terms of the instantaneous abundances of the neu- 
tral species. This simplification is based on the 
fact that the "ion chemistry" timescale is much 
shorter than the depletion timescale (see Caselli 
et al. 2002b and Crapsi et al. 2004for details). 
Here we included the multiply deuterated forms 
of , without any distinction between the ortho 
and para forms of molecular hydrogen and H 2 D + 
and without taking into account the so-called 
back-reactions between, e.g., ortho-H2 and ortho- 
H 2 D+ (Gcrlich, Herbst, & Roucff 2002), which 
limit the deuterium fractionation, as discussed in 
Walmsley, Flower, & Pineau des Forets (2004). 
However, we adopted the new value of the rate 
coefficient for the proton-deuteron exchange reac- 
tion H++ HD -> H 2 D++ H 2 (3.5xl0~ 10 cm 3 s" 1 ; 
Gerlich, Herbst, & Roueff 2002). The other pa- 
rameters used are: the cosmic-ray ionization rate 
C = 1.3xl0~ 17 s" 1 ; the CO, N 2 , and O binding 
energies: E(CO) = 1210 K, E(N 2 ) = 800 K, and 
E(0) = 650 K 3 respectively; and the lower cutoff 
radius of the MRN (Mathis, Rumpl, & Nordsieck 
1977) distribution a m i n = 5xl0~ 6 cm. 

The different model predictions, denoted as m, 
n 2 , ri3, n 4 and n 5 in Figure 5, refer to model 
spherical clouds with different density structures, 
in (rough) analogy with the density structure of 
the model cloud undergoing infall, described by 
Ciolek & Basu (2000, hereafter CB00), at succes- 
sive evolutionary stages (from t\ = 2.27 Myr to 
t 5 = 2.684 Myr: see CB00). We remark that our 
calculations assume the steady state and thus our 
model makes use of the density profile only, and 
not of the time scale. We also note that the CB00 
model has cylindrical (rather than spherical) sym- 
metry and we have neglected this fact assuming 
the same dependence on spherical radius as that 



All the adopted parameters are the same as in Crapsi et 
al. (2004), with the exception of E(0), which is now 100 K 
lower than before. The decrease in this parameter is neces- 
sary to keep the deuterium fractionation at a similar level 
as in Crapsi et al. (2004), for the particular case of L1521F, 
after the inclusion of all the multiply deuterated forms of 
Hjj~. See also Caselli et al. (2002b) for the importance of 
atomic oxygen in this simple chemical network. 
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on cylindrical radius in CBOO. We assume that the 
density profiles are given by the analytical formula 
of Tafalla et al. (2002), with different values of the 
parameters (namely, the central density, the ra- 
dius of the inner flat region, and the asymptotic 
power index) to approximately reproduce the den- 
sity profiles of Figure la of CBOO. In particular, 
the central densities of the cores at each time are 
given by 10 j x 4.37 x 10 3 cm~ 3 (for j = 1, 5) as 
in CBOO. The abundance profiles obtained in this 
way have then been converted into column densi- 
ties via integration along the line of sight and suc- 
cessive convolution with a two-dimensional Gaus- 
sian, simulating observations with an HPBW of 
22". From Figure 5 we note that the deuterium 
fractionation in N2H + and the observed CO deple- 
tion factor (/d) are predicted to increase with core 
evolution and actually explain the observed trend, 
although some scatter is expected in the data, 
given the different environments where they are 
immersed and their different formation histories. 
In particular, we expect magnetic field strength, 
amount of turbulence, external radiation field and 
external pressure to influence the evolution of a 
core. Hence two cores with the same mass and 
same age could have different degrees of evolution 
depending on the characteristics of the surround- 
ing star-forming region. We note that most of the 
"pre-stellar" cores are found close to the n 3 re- 
sult, which, in the CBOO model, is the model of 
best agreement with the measured density profile 
of L1544. Therefore, we can conclude that the 
"pre-stellar core" condition, 7V(N 2 D+)/7V(N 2 H+) 
> 0.1, is consistent with values of central densities 
n (H 2 ) > 10 6 cm~ 3 , in agreement with H 2 central 
densities measured from the dust in Section 4.3 for 
these cores. 

We also searched for correlations between the 
depletion factor and the H 2 column density and 
found a positive correlation (87% correlation coef- 
ficient; 90% for Taurus only), confirming the gen- 
eral chemical scenario about CO depletion (see 
Dalgarno & Lcpp 1984; Roberts & Millar 2000; 
Bacmann et al. 2002). Several authors have inves- 
tigated the chemical structure of pre-stellar cores, 
coupling together dynamical and chemical evolu- 
tion with detailed models (Bergin & Langer 1997; 
Aikawa et al. 2001; Li, et al 2002; Shematovich, 
et al 2003). Deuterium and singly deuterated 
species were, however, only included in Aikawa et 



al. (2001), who were able to reasonably reproduce 
the observed molecular D/H ratios toward L1544, 
although the column density of N 2 H + was under- 
estimated. Aikawa, Ohashi, & Herbst (2003) over- 
came this problem by adding surface chemistry to 
their models. Surface processes significantly en- 
hanced the production of the parent molecule N 2 , 
increasing the N 2 H + column density to the ob- 
served values. 

More recently, Roberts, Herbst & Millar (2004) 
studied the chemistry of pre-stellar cores includ- 
ing multiply deuterated species but neglecting the 
dynamical evolution. For the particular cases of 
L1544 and OphD they found 7V(N 2 D+)/7V(N 2 H+) 
column density ratios significantly (factor of about 
5) larger than the observed values, suggesting that 
the deuterium fractionation process is probably 
too efficient in their models. Possible causes of 
this disagreement are: (i) the exclusion of the so- 
called back reactions of deuterated isotopologues 
of H^ with ortho-H 2 (Gerlich, Herbst, & Roueff 
2002) which may lower the molecular D/H ratios 
(Walmsley, Flower, & Pineau des Forets 2004); (ii) 
the use of large (radius of 0.1 /im) dust grains in 
the chemical network, which underestimates the 
recombination of molecular ions compared to mod- 
els where a population of smaller grains is present 
or where an MRN size distribution of dust par- 
ticles is considered; (iii) the use of UMIST rates, 
which produce Hg and analogues more efficiently 
than the "New Standard Model" from Ohio State 
University (see Roberts, Herbst & Millar 2004, 
for details). 

5.1.2. Deuterium fractionation vs. _ff 2 central 
density, dust equivalent radius, N 2 H + line 
width and N 2 H + line skewness 

In theoretical simulations, we expect an in- 
crease of the central n(H 2 ) density and a decrease 
of r7o, the equivalent radius of the flat region, with 
ongoing contraction (e.g., Lizano & Shu 1989; Fos- 
ter & Chevalier 1993; Ciolek & Mouschovias 1995; 
Li 1999). Those two parameters can be consid- 
ered as indicators of "dynamical evolution" ; here 
we cross-correlate them with the deuterium frac- 
tionation. 

From Figure 6 panel a) one can see that there is 
a positive correlation between 7V(N 2 D + )/A^(N 2 H + ) 
and n(H 2 ) apart from OphD and L1521F which 
show opposite behaviours from the rest of the 
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cores. Six of our pre-stellar cores candidates 
have H2 volume densities above the average 
(5.1 x 10 5 cm -3 ), and all of these show increases in 
deuterium fractionation above the average. Our 
estimate of the volume density depends on the as- 
sumed spherical geometry: hence, we also consid- 
ered the correlation between the deuterium frac- 
tionation and the H 2 column density. While the 
general trend remained the same, the peculiarity 
of OphD and L1521F becomes more marked. 

We check here the relative variation of the flat- 
tened region in the 13 cores for which we have 
a 1.2-mm map. Having treated all the cores in 
a consistent way allows us to consider the differ- 
ences in r7o real and not due to differences in the 
adopted fitting model. 

The correlation between deuterium fractionation 
and is shown in Figure 6 panel b). A close 
relationship of the flattened radius with the deu- 
terium fractionation is not found. 
OphD shows the biggest variation from this pic- 
ture, but one should notice the presence of a sec- 
ond peak inside its 70% contour. It is thus possible 
that the OphD equivalent radius is overestimated. 
A similar argument holds for L183. Limiting the 
sample to the Taurus cores significantly tightens 
the correlation, increasing the correlation coeffi- 
cient from 1% to 64%. Once again, this result 
suggests that cores embedded in different environ- 
ments (e.g. external radiation field) hardly rep- 
resent a homogeneous sample in which to study 
evolutionary trends. 

Caselli et al. (2002a) showed that in L1544 there 
was a definite trend for the N2H 4 " and N2D 4 " line 
widths to increase towards the high-density region. 
This result was interpreted as a sign of increasing 
infall activity towards the center. The same trend, 
although less marked, is seen in L1521F (Crapsi et 
al. 2004). 

In this fashion, we take here the N2H" 1 " line width 
as an indicator of infall activity and thus prox- 
imity to the critical state for dynamical collapse 
and we check if the cores that had larger line 
widths arc the same as those with the largest 
iY(N 2 D+)/7V(N 2 H+) ratio, to search for links be- 
tween physical and chemical evolution. 
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Fig. 6. — Deuterium fractionation correlations. In all panels: dotted lines are the arithmetic means of 
the distributions (except in panel d)); cores belonging to the Taurus Molecular Cloud are marked by a 
square; pre-stellar cores candidates are circled; the values were evaluated at N 2 D + peak position. Panel a): 
Deuterium fractionation vs. H 2 central volume density. The n(H 2 ) values are taken from the literature (see 
Section 4.3), with the exception of TMC2, L492 and B68, estimated in this work. Panel b): Deuterium 
fractionation vs. equivalent radius of the flattened region tjq is defined as \J aream/it, with area^ as 
the area within the 70% contour in the 1.2-mm map. Panel c): Deuterium fractionation vs. N 2 H+(l-0) 
line width evaluated from the simultaneous fit of all the 7 hyperfine components. Panel d): Deuterium 
fractionation vs. N 2 H + (l-0) line asymmetry. The vertical dotted line separates the red-skewed spectra from 
the blue-skewed ones (the arithmetic mean is 0.07). 
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The result is in Figure 6 panel c). In our sample 
we do not see any correlation for larger line width 
in the nucleus of cores with larger deuterium frac- 
tionation. In particular, we note that L328 shows 
broad lines but low N (N 2 D+)/7V(N 2 H+) and that 
OphD has the opposite behaviour. The two broad- 
est line cores, L328 and L429, have spectra with 
both flat tops and blue-skewed asymmetry (see 
Sect. 4.5), suggesting that their lines could be 
broadened for reasons other than infall. We at- 
tempted to avoid the environmental differences, 
as a source of scatter for the correlation, limiting 
the sample to the Taurus cores only; in this case 
the correlation coefficient increases from 22% to 
40%. 

We point out that the above correlation does not 
make any distinction between cores that do and 
do not show evidence of central infall (through the 
N2H + (l-0) line width broadening toward the dust 
peak position, or with the presence of "red shoul- 
ders"), so that the different line widths observed 
in the cores of the present sample may simply 
arise from different amounts of turbulence in their 
interiors or other types of motions (e.g., expansion 
or oscillations) along the line of sight. 

We also try to correlate the line skewness 
with the iV(N 2 D+)/7V(N 2 H+) ratio in Figure 6 
panel d). Also in this case there is a general 
agreement with the proposed idea that cores with 
higher deuterium fractionation also show signs 
of kinematical evolution. The main exceptions 
in this case are given by L429 on the one hand 
and TMC1C and L492 on the other. It is in- 
teresting to note that the correlation between 
7V(N 2 D + )/7V(N 2 H + ) and the line asymmetry is 
more tight than that with the line width. 

5.1.3. Deuterium fractionation and extended in- 
fall 

Lee, Myers, & Tafalla (1999, 2001) and Lee, 
Myers, & Plume (2004) in a series of papers con- 
ducted a search for infall motions towards a sam- 
ple of 70 starless cores. Infall motions were identi- 
fied through the presence of velocity shifts between 
high density tracers (N 2 H + (l-0) probing densities 
up to a few times 10 6 cm -3 ) and low density trac- 
ers (CS(2-1) tracing the gas up to ~ 10 4 cm~ 3 ) 
and through the observation of double-peaked self- 
absorbed optically thick lines whose blue peak was 



brighter than the red peak (Lee, Myers, & Plume 
2004). Given the CS depletion in the inner core, 
those measurements gauge the infall of the exter- 
nal envelope only. 

Twenty-four cores in our sample were observed 
also in the Lee ct al. papers (2001; 2004), in- 
cluding our chemically evolved candidates L1544, 
L1521F, L694-2, L183, L429, TMC2 and OphD 
(there known as L1696A). Considering now only 
this common subsample, Lee, Myers, & Plume 
(2004) found that 12 of these cores show in- 
fall signatures in at least three of the seven 
tests they performed on CS, DCO+ and N 2 H+ 
lines, thus they are considered strong infall can- 
didates. All our chemically evolved candidates 
(iV(N 2 D+)/iV(N 2 H+)> 0.1) (except OphD, which 
was not studied with all the techniques by Lee, 
Myers, & Plume 2004) were identified as strong in- 
fall candidates by the Lee, Myers, & Plume (2004) 
analysis. This seems to show that our chemically 
evolved candidates, the majority of which show 
clear signs of central infall, are a subsample of the 
cores undergoing extended contraction. 

We also performed a direct comparison between 
one of the infall indicators in Lee, Myers, & Plume 
(2004), <5V DCO + {= (V DCO+ - V N2H+ )/AV N2 H+), 
and our N 2 H + (l-0) skewness on a common sub- 
sample of 15 cores. We find a very good agree- 
ment between those two parameters with only two 
cases in which they give a different indication: 
L429 which has significant outflow signature in 
the N 2 H + (l-0) skewness but marginal infall in- 
dication in <5Vrjco+7 an d L183, which has a very 
strong outflow behaviour in <5V DCO + whereas it 
shows an infall profile in our N 2 H + (1- 0) spectrum. 
For the 13 remaining cores the two indicators cor- 
relate with a 70% correlation coefficient. 

5.2. Dynamical parameters correlations 

We cross-correlated our "kinematic activity" in- 
dicators, i.e., the N 2 H + (l-0) line width, the ra- 
dius of the flattened region, the H 2 column den- 
sity, and the N 2 H + (l-0) line asymmetry, to see 
whether they give consistent information. 
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Fig. 7. — Correlations between "dynamical evolution" indicators. In all panels: dotted lines are the arith- 
metic means of the distributions (but in panel c)); cores belonging to the Taurus Molecular Cloud are marked 
by a square; pre-stellar cores candidates are circled; the values were evaluated at N 2 D + peak position. Panel 
a): N 2 H + (l-0) line width vs. equivalent radius of the flattened region is defined as y^areaia/n, with 
arearo as the area within the 70% contour in the 1.2-mm map. Panel b): N2H + (l-0) line width vs. H2 vol- 
ume density. The n(H 2 ) values are taken from the literature (see Section 4.3), with the exception of TMC2, 
L492 and B68, estimated in this work. Panel c): N 2 H + (l-0) line width vs. N 2 H + (l-0) line asymmetry. 
The vertical dotted line here separates the red-skewed spectra from the blue-skewed ones (the arithmetic 
mean is 0.07). Panel d): N 2 H + (l-0) line width vs. aspect ratio of the 50% contour of the 1.2-mm emission. 
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The result is encouraging. We can in fact see in 
panel a) of Figure 7 that rra is anti-correlated 
with N2H + (l-0) line width (correlation coefficient 
67%; 86% limiting to Taurus cores); note that 
even cores that in the other plots do not behave 
like the rest of the sample (OphD and L328) lie 
here in the common trend. 

In panel b) of Figure 7 we find that, with the no- 
table exceptions of L328 and L429, the N 2 H+(l-0) 
line width correlates fairly well with the H 2 cen- 
tral density (correlation coefficient 20% for the to- 
tal sample; 64% without L328 and L429; 94% for 
only the Taurus cores). 

In the picture given in Section 5.1.2, the cores 
with strong motions along the line of sight should 
show stronger asymmetries and broader lines to- 
wards the core peak. In Figure 7 panel c) we plot- 
ted the relation between the width and the skew- 
ness of the N2H 4 " line. Although again a tight cor- 
relation is not found, we can see that there is a gen- 
eral trend of increased skewness in broader spec- 
tra. Remarkable exceptions are constituted by 
OphD and TMC1-C, having narrow lines strongly 
skewed towards the red, and by L429 and L328, 
which have the broadest lines in the sample but 
do not show any infall asymmetry. Also in this 
case, limiting the sample to Taurus cores, the cor- 
relation coefficient increases from 2% to 21%. 

5.2.1. Line width and aspect ratio 

The fact that the majority of the cores in our 
sample (sec Figure 14 in Section B) show large 
departures from spherical symmetry has consid- 
erable significance for our understanding of the 
dynamics of these objects. Core shape could for 
example reflect the magnetic field structure and 
indeed several of the cores in our sample show lin- 
ear polarization of around 5%-10% (e.g., Crutcher, 
Nutter, Ward-Thompson, & Kirk 2004), suggest- 
ing magnetic fields of order 100 fiG. However, the 
observed structures could also be a consequence 
of dynamical interactions with the core surround- 
ings as is suggested by the observed "cometary" 
shapes. It is natural therefore to search for kine- 
matic evidence in our N 2 H+ maps that either ex- 
cludes or favors one of these scenarios. 

There are various possible kinematic probes, in- 
cluding the observed velocity field, but the sim- 



plest is to check for a correlation between line 
width and aspect ratio. Here the interpretation 
depends considerably on whether the true core 
shape is "close to oblate" or "close to prolate" (see, 
e.g., Jones & Basu 2002; Curry 2002, who have 
attempted to decide this on the basis of the statis- 
tics of observed aspect ratios) . Present results are 
inconclusive and we conclude that one must con- 
sider both possibilities. 

In either case, if cores are magnetically domi- 
nated, motion along the field lines will be facili- 
tated. We would expect that in such cases, mo- 
tions will predominantly be parallel to the axis of 
symmetry and hence for oblate cores, one might 
expect the largest line widths for face-on circu- 
lar sources and for prolate cores for end-on (also 
circular) objects. On the other hand, ambipolar 
diffusion necessarily occurs across field lines and 
this will have the opposite effect. Thus, for exam- 
ple, edge-on oblate cores will have the largest line 
width. In fact, CB00 have provided predictions 
for how the velocity field due to the latter process 
might appear in L1544 and Caselli et al. (2002a) 
showed that the observations were roughly consis- 
tent with the model for this object. 

However, this does not seem to be the case 
in L1521F (Crapsi et al. 2004) and also not in a 
more general sample of cores as shown in Figure 7 
panel d), where we plot the N2H + (l-0) line width 
against aspect ratio from the millimeter contin- 
uum maps. There is no clear trend in this plot, and 
clearly a larger sample is needed. However, our re- 
sults do not give evidence that in general there are 
increased line widths towards elongated core nu- 
clei as one might expect on the basis of models 
such as those of CB00. On the other hand, our re- 
sults also do not seem to favor the models of Curry 
& Stahlcr (2001) which suggest that initially elon- 
gated structures should become more spherical as 
they evolve owing to ambipolar diffusion. We con- 
clude that it would be useful to have more detailed 
predictions for model pre-stellar cores of the evo- 
lution with time of both the line width and the 
velocity field. 

5.3. Associated properties 

We have shown in this section a series of prop- 
erties that we think are associated with an evolved 
stage of starless cores either because of enhanced 
chemical activity or because of an advanced kinc- 
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Table 9 

Chemical and kinematical evolution "probes" in the 13 cores with full data set (N 2 H + , 

N 2 D + , C 18 AND 1.2-MM EMISSION). 



core 


LI 544 


L15211 


L694-2 


L429 


L183 


OphD 


TMC2 


L492 


L328 


L1495 


L15171 


00 

m 


L1498 


Af(N 2 D+) > 1.0 x 10 12 cm" 2 


1 


1 


1 


1 


1 


1 


1 




















A r (N 2 H+) > 8.5 x 10 12 cm" 2 
lN2D+] - 10 

[N 2 H+] > U - iU 


1 


1 


1 


1 


1 





1 


1 


1 


1 


1 


i 





1 


1 


1 


1 


1 


1 


1 




















f D (CO) > 10.2 


1 


1 


1 


1 


1 


1 


1 




















n(H 2 ) > 5.1 X 10 5 cm- 3 


1 


1 


1 


1 


1 


1 























a Vn 2 h+ > 0-25 km s- 1 


1 


1 


1 


1 











1 


1 














infall asym. (skewness > 0) 


1 


1 


1 





1 


1 





1 





1 











r 70 < 4800 AU 


1 


1 





1 


1 











1 


1 


1 


1 





total 


8 


8 


7 


7 


7 


5 


4 


3 


3 


3 


2 


2 






Note. — For each value, but the skewness, the thresholds are given by the arithmetic mean 
value for the sample. 



matical process. 

Even though we did not find tight one-to-one 
correlations between the numerical values it ap- 
pears that these properties globally isolate a group 
of cores that look peculiar with respect to the rest 
of the sample. We summarize these characteristics 
of the evolved starless cores, or pre-stellar cores, 
in Table 9. L1544 and L1521F were found to be 
more evolved than the "average" core in each test 
we performed, while L183, L429 and L694-2 failed 
just one test and can be considered strong can- 
didates for pre-stellar cores as well. Other cases, 
such as OphD and TMC2, are more doubtful. 

In the presence of an unbiased survey and as- 
suming that all the cores will eventually form a 
star, one could attempt to estimate the duration 
for the pre-stellar core phase by dividing the num- 
ber of pre-stellar core candidates by the total 
number of objects in the sample. Since our main 
selection criterion stated that only the sources 
with stronger N 2 H + (l-0) observations (in the Lee, 
Myers, & Tafalla 2001samplc) or stronger dust 
emission-absorption (in the Bacmann et al. 2000 
sample) were observed, we cannot consider our 
survey unbiased. Nevertheless, taking the num- 
ber of N2H + -detected cores within 250 pc (49) in 
the Lee, Myers, & Tafalla (1999) sample as repre- 
sentative of the total number of the N 2 H + -bright 
cores and supposing that L1544, L1521F, L183, 



L429, L694-2, OphD and TMC2 are the only pre- 
stcllar cores observable in this sample, we can de- 
rive a duration of the pre-stellar phase equal to 
7/49, or —14% of the lifetime of a core since it 
started to be observable in N 2 H + . Taking the 
smallest central density observed in our sample 
(B68: n(H 2 ) ~ 8 x 10 4 enr 3 ) as the threshold 
for N 2 H + emission, we can derive a lifetime of an 
N 2 H+ -bright core of —0.4 Myr (using the CB00 
ambipolar diffusion models); hence, the pre-stellar 
core phase should last around 60,000 yr. An alter- 
native estimate can be made evaluating the free- 
fall time for the threshold density for a pre-stellar 
core (n(H 2 ) — 5.1 x 10 5 ) resulting in a duration of 
-60,000 yr. 

The pre-stellar core lifetime can also be esti- 
mated from the comparison with the estimated 
time spent as a Class protostar (~ 2 x 10 5 yr 
Visser, Richer, & Chandler 2002). From Andre, 
Ward-Thompson, & Barsony (2000) and Visser, 
Richer, & Chandler (2002), we derive a total num- 
ber of 11 Class protostars within 250 pc. This, 
compared to the seven evolved cores found in this 
work, gives a duration of the pre-stellar phase 
of —120,000 yr, longer than the estimates above. 
However, we feel that this discrepancy merely re- 
flects the great uncertainties in both estimates. 
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6. Conclusions 

We collected a set of N 2 H+, N 2 D+, C 18 
and 1.2-mm continuum observations using the 
IRAM 30-m for a sample of 31 starless cores, com- 
pleting literature data with new observations. We 
retrieved from these data chemical and kinemat- 
ical parameters supposed to discriminate young 
from more evolved cores such as: deuterium frac- 
tionation, integrated CO depletion factor, H 2 den- 
sity, N 2 H+ line width and line asymmetry. Our 
conclusions are summarized as follows: 

1. We attempted several correlations between 
the chemical evolution indicators (deuterium frac- 
tionation, integrated CO depletion factor and 
N 2 H + and N 2 D + column density) and properties 
of dynamically evolved cores (such as H 2 volume 
density, line width and asymmetry and size of the 
flattened region). Although we did not find very 
tight dependencies between these quantities, we 
recognized that all these properties, as a whole 
more than taken singularly, draw from the total 
sample those cores with an enhanced state of evo- 
lution. As can be seen from Table 9, we find that, 
compared to the total sample, L1544, L1521F, 
L183, L429, L694-2, OphD and TMC2 globally 
show: 

• higher N 2 H + column density 

• higher N 2 D + column density 

• higher iV(N 2 D+)/A(N 2 H+) abundance ra- 
tio 

• higher H 2 column density 

• higher integrated CO depletion factor 

• broader N 2 H+ lines 

• smaller flattened radius 

• stronger infall asymmetries 

with some exceptions in each category. 

2. Limiting the sample to the cores belonging 
to the Taurus Molecular Cloud tightens the cor- 
relations, probably by diminishing possible effects 
due to environmental differences such as magnetic 
field strength, amount of turbulence, external ra- 
diation field and external pressure. 

3. The picture emerging from the observations 
is consistent with predictions of chemical models 



accounting for molecular freezeout, where cores 
with higher central densities show a lower frac- 
tion of CO in the gas phase and a correspondingly 
larger deuterium enrichment in the core nucleus. 

4. N 2 D+ and N 2 H+ emissions, unlike C 18 0, 
peak very close to the dust continuum max- 
ima in the mapped cores, confirming the fact 
that those species probe better the conditions 
of the high density gas. Moreover, the maps of 
N 2 D + and N 2 H + (3-2) were found to be systemat- 
ically smaller than the N 2 H + (l-0) ones in all the 
mapped cores, consistent with the idea that these 
trace better the high-density nucleus. We note 
that central H 2 densities evaluated from the dust 
emission were found systematically larger than 
those calculated from LVG modeling of N 2 H + and 
N 2 D + data, suggesting that depletion of N 2 could 
possibly occur at the highest density peak, al- 
though finer resolution data are needed to confirm 
this effect. 

5. The decrease of line-width towards the 
edge observed in L1544 (Caselli et al. 2002a) and 
L1521F (Crapsi et al. 2004) has been found in only 
one other core (L694-2) with supposedly the same 
degree of evolution. 
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A. Maps 



Figures 8 to 13 show maps towards L183, OphD, L429, TMC2, L492 and L694-2, which, together with 
L1544 and L1521F, stand out from the rest of the sample as having very bright N 2 D+(2-l) lines. Similar 
maps for L1544 can be found in Caselli et al. (2002a) and Tafalla et al. (2002) and for L1521F in Crapsi et 
al. (2004). The N2H" 1 " and N2D" 1 " data shown in Figures 8-11 were all convolved with a 26" beam in order to 
reach the same spatial resolution as the N 2 H + (l-0) data and make the comparison easier. The N 2 H + (l-0) 
map of TMC2 (Figure 13) was kept with the original FCRAO resolution (54"), while we have not smoothed 
the N 2 D + (2-l) map of L492 (Figure 12) because of the small extension of it. Similarly, all the 1.2-mm maps 
were smoothed to a 22" resolution, the beamsize of the IRAM 30-m at the C 18 O(l-0) frequencies. 

The most evident feature of Figures 8 to 13 is that the C 18 does not trace the dust continuum while the 
N 2 H + and N 2 D + emission morphology is similar to that of the dust, as already shown by several authors 
(e.g. Caselli et al. 1999; Bergin et al. 2001; Tafalla et al. 2002) . While this is true for L1544, L1521F, 
L429, L694-2 and OphD, there are differences between the 1.2-mm map and the N 2 H + -N 2 D + maps for L183 
and TMC2. For example, the L183 dust emission peaks 30" to the south of N 2 H + (l-0) and, especially, 
N 2 D+(2-l). 

TMC2 shows an even more complex situation; in fact, the 1.2-mm maps revealed much more structure 
than the N 2 H + (l-0) map of Caselli, Benson, Myers & Tafalla (2002) led us to believe. Whereas we found a 
local 1.2-mm peak next to the N 2 H + (l-0) one, we also found brighter peaks of the dust emission in regions 
where the N 2 H+(l-0) is fainter. Note that the N 2 H+ map in Figure 13 has a lower angular resolution 
(HPBW = 54") than in the other sources; smoothing the dust continuum map to this resolution does not 
help to reconcile the two emissions, so we present here the 1.2-mm map with the same resolution as in the 
other cores (i.e. 22"). We also add in Figure 13 the position of a few radio stars and of a couple of IRAS 
point sources; the presence of these objects could explain the discrepancy between the N 2 H + and the dust 
maps. 

Another common feature clearly seen throughout Figures 8 to 13 is that the N 2 D + maps are systematically 
smaller than the N 2 H + ones (see Table 10). This could indicate that the deuterated species are better tracers 
of the core nucleus. Similar behaviour is seen when the higher transition maps are compared with those at 
lower J; this may be caused by the increase in density towards map center. 

In OphD (see Figure 10), we note the presence of two nuclei embedded in the same envelope; both have 
a similar brightness in the 1.2-mm continuum although they have a very different intensity in the N 2 H + (1- 
0), N 2 D+(2-l) and N 2 D+(3-2) (the N 2 H+(3-2) has poor S/N towards the northeast core). The relative 
intensity of the two nuclei decreases going from the lowest density tracer (N 2 H+(l-0)) to the highest density 
one (N 2 D + (3-2)), suggesting that the southwest core is more centrally peaked than the northeast one. 



Table 10 

n 2 h+, n 2 d+ and 1.2-mm maps full width at half maximum. 



core 


1 . 2mm 


N 2 H+(l-0) 


N 2 H+(3-2) 


N 2 D+(2-l) 


N 2 D+(3-2) 


L183 


185" x 67" 


217" x 92" 




93" x 72" 




OphD 


129" x 70" 


170" x 70" 




93" x 62" 


67" x 65" 


L492 


180" x 100" 


98" x 78" 








L429 


102" x 53" 


91" x 78" 


72"x53" 


60" x 56" 


43" x 41" 


L694-2 


77" x 55" 


85" x 67" 




68" x 60" 
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Fig. 8.— L694-2 maps in 1.2-mm continuum, C 18 O(l-0), N 2 H+(l-0) and N 2 D+(2-l) . N 2 H+ and N 2 D+ 
data were smoothed to a common resolution of 26" and rcgridded; the reference position is the one reported 
in Table 1. The 1.2-mm map (M. Tafalla et al. 2005, in preparation) was smoothed to a resolution of 22". 
Contour levels are 30% to 90% by 15% of the peak for each map. Peak values are 185 mJy/22" beam for 
the 1.2-mm map and 5.16, 1.26 and 1.76 K km s" 1 for N 2 H+(l-0), N 2 D+(2-l) and C 18 O(l-0) respectively. 
Available 40"x40" maps of N 2 H + and N 2 D + (3-2) are not reported because of poor S/N outside the peak. 



Fig. 9.— L183 maps in 1.2-mm continuum, C 18 O(l-0), N 2 H+(l-0) and N 2 D+(2-l). Same as in Figure 8 
but for the peak values: here 239 mJy/22" beam for the 1.2-mm map and 4.31, 1.65, and 2.5 K km s _1 for 
N 2 H + (1 -0), N 2 D + (2-l) and C 18 O(l-0) respectively. The 1.2-mm map has been taken from Pagani et al. 
(2003). Available 30"x60" maps of N 2 H+and N 2 D+(3-2) arc not reported because of poor S/N outside the 
peak. 



Fig. 10.— OphD maps in 1.2-mm continuum, N 2 H+(l-0) and (3-2) and N 2 D+(2-l) and (3-2). Same as in 
Figure 8 but for the peak values: here 195 mJy/22" beam for the 1.2-mm map and 3.81, 0.49, 1.75 and 0.60 
K km s- 1 for N 2 H+(l-0) and (3-2), N 2 D+(2-l) and (3-2) respectively. The 1.2-mm map has been taken 
from Ward- Thompson, Motte & Andre (1999). The RMS of the N 2 H+(3-2) increases significantly towards 
north east, thus contours above -20" have to be taken with caution. 



Fig. 11 — L429 maps in 1.2-mm continuum, N 2 H+(l-0) and (3-2) and N 2 D+(2-l) and (3-2). Same as in 
Figure 8 but for the peak values: here 210 mJy/22" beam for the 1.2-mm map and 6.12, 0.63, 2.07and 0.69 
K km s- 1 for N 2 H+(l-0) and (3-2), N 2 D+(2-l) and (3-2) respectively. The 1.2-mm map has been taken 
from Bacmann et al. (2000). 



Fig. 12.— L492 maps in 1.2-mm continuum, C 18 O(l-0), N 2 H+(l-0) and N 2 D+(2-l). Same as in Figure 8 
but for the peak values: here 114 mJy/22" beam for the 1.2-mm map and 3.96, 0.31 and 1.69 K km s _1 for 
N 2 H+(l-0), N 2 D+(2-l) and C 18 O(l-0) respectively. 



Fig. 13. — TMC2 maps in 1.2-mm emission smoothed to 22" (contours in both panels) and C 18 O(l-0) (gray 
scale in the top panel) from new 30m observations and the N 2 H + (l-0) map (gray scale in the bottom panel) 
taken by Caselli, Benson, Myers & Tafalla (2002) with FCRAO (HPBW 54"; reference position here is the 
same as in Table 1). Contour levels starts from 30% of peak and increase by 15% of it. Peak values are: 
138 mJy/22" beam for the 1.2-mm map, 2.12 K km s" 1 for the C 18 observations. The N 2 H+(l-0) intensity 
peak is 3.6 K km s _1 in Tmb temperature. In the bottom panel the positions of seven young stars (in order 
of ascending RA: JH93, JH90, JH91, JH92, JH94 - Wendker 1995 ; FY Tau, FZ Tau - Nurmanova 1983) arc 
denoted by stars. Two IRAS point sources visible in the same field are denoted by a filled square. 
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B. Core shapes and kinematics 

We underline a marked tendency of the starless cores for having the dust emission peak shifted relative 
to the centroid of the cores, suggesting that the cores either have been compressed by external agents or 
have moved relatively to the surrounding environment (but see Walsh, Myers & Burton 2004). In Figure 14 
we attempt a classification of the shape of the starless cores evaluating a parameter able to quantify the 
degree of "cometariness" or how much the dust peak is shifted relative to the lower contours of the map. We 
evaluated the centroid of the map as the mean position of the points within the 50% contour and defined 
the degree of "cometariness" as the distance between this and the dust peak in units of r-jo, the equivalent 
radius of the core. OphD and TMC2 were not included in this classification since they harbour multiple 
peaks within their 50% contour (see Figures 10 and 13); thus, the "cometariness" parameter is meaningless. 
From this small subsample we do not notice any clear correspondence between a more marked cometary 
shape and other kinematical properties of the core (e.g. H 2 density, N 2 H + (l-0) line-width, or r-jo). 



Fig. 14. — The 1.2-mm continuum maps of 12 cores in our sample. Light contours are spaced from 10% to 
90% by 10% of the peak, while the 98% and 50% contours are represented by a thick line. The cross marks 
the mean position of the points within the 50% contour. The cores are ranked for increasing "cometariness" . 
This parameter is reported at the top right corner of each map and is calculated as the distance between the 
mean point and peak in units of rm . TMC2 and OphD were not included in this classification since they 
have multiple peaks within their 50% contour, so this technique failed for them. 



Fig. 15.— N 2 H+(l-0) velocity gradients towards the cores L694-2, OphD, L183, L429 and L492. The 
N 2 H + (l-0) data were regridded in a regular pattern. The integrated intensity maps are shown in greyscale 
(see Figures 8 - 12). Local velocity gradients in the adjacent nine points are represented by an arrow pointing 
in the direction of increasing velocity and with the length proportional to the magnitude (1 km s _1 pc -1 
corresponds to a 6"-long arrow). Note that the angular scale of each map is the same. 



Using the N 2 H + (l-0) velocities from the hyperfine structure fitting, we found strong velocity variations in 
L183 (0.23 km s _1 in a 0.16 pc length) and OphD (0.19 km s _1 in a 0.12 pc length), while a smaller variation 
(0.07 km s _1 in a 0.1 pc length) was found in L694-2 and constant velocity fields were seen in the inner 
0.08 pc of L429 and L492. We evaluated the total velocity gradients of these cores under the assumption of 
solid body rotation (see Goodman, Benson, Fuller & Myers 1993). The magnitudes and the direction of the 
velocity gradients are reported in columns (1) and (2) of Table 11. 

These calculations agree with the velocity variations seen across the cores. 
Local velocity gradients, evaluated in 3 x 3 adjacent points, were also calculated towards these cores (see 
Caselli et al. 2002a, for details on the procedure). The analysis was performed on regridded data and 
results are shown in Figure 15 and Table 11. L183 and OphD present strong motions and a regular velocity 
field. L429 and L492 show very little motion in the inner parts but larger gradients at the edge of the map 
with arrows pointing in different directions on different sides of the cores, explaining the absence of velocity 
variation in the inner part and the small total velocity gradient. Finally, L694-2 shows the smallest local 
gradients with a regular velocity field only in the inner part. 

Using the N 2 H + (l-0) and N 2 D + (2-l) maps, we also checked for the presence of line width gradients in 
L694-2, OphD, L183, L429 and L492 in the same fashion as for L1544 (Caselli ct al. 2002a), B68 (Lada, 
Bergin, Alves, & Huard 2003) and L1521F (Crapsi et al. 2004). From figure 16 we find that only L694-2 
shows the line width decrease that was evidenced in L1521F and L1544, while a marginal similarity is found 
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Fig. 16.— N 2 H+(l-0) intrinsic line width gradients within the cores L694-2, OphD, L183, L429 and L492. 
The angular offset was evaluated from the N 2 D + (2-l) peak position (see Table 3). 
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in the first 50" around OphD and in the northern part of L492. The opposite behaviour is seen in L429 and 
L183 (as found in B68). The N2D + (2-l) line widths evidenced the same behaviour as in figure 16. 

C. Density determination in the LVG approximation 

The Large Velocity Gradient (LVG) approximation simplifies the radiative transfer treatment assuming 
that a photon produced in a contracting or expanding cloud is reabsorbed locally and that the local velocity 
gradient is much larger than the ratio of line width to cloud size. In addition, one normally assumes 
homogeneous and isothermal conditions. Assuming a temperatutre of 10 K, we have used this approximation, 
with collisional rates from Green (1975) to compute level populations for N 2 H+ and N 2 D + and to predict 
line intensities for these species as a function of density and species column density. We then compared the 
model predictions with the observed N 2 H + (3-2)/N 2 H + (l-0) intensity ratio and the N 2 H + (l-0) integrated 
intensity (or the N 2 D+(3 2)/N 2 D+(2 1) and N 2 D+(2 1) intensity) exploiting the fact that the former is 
mainly sensitive to n(H 2 ) and the latter to the N 2 H + (N 2 D + ) column density. Our results are extremely 
sensitive to the intensity of the two 3-2 transitions and we hence confine our discussion to the eight cores 
with S/N> 5 in the 1 mm lines. An extra constraint comes from the total r of N 2 H + (l-0) (N 2 D + (2-l)). 
Note that, given the small opacities found in N 2 D + (2-l) (and in many N 2 H+(l-0) spectra), r has little 
influence on the n(H 2 ) determination whereas it is a useful constraint for column density. 

In Table 12 we summarize the results of the LVG calculations. The ratio of the N 2 H+to N 2 D+ column 
densities inferred from LVG and the relative value found in the CTEX approximation have an arithmetic 
mean value of 0.9 and a dispersion of 25%. Thus, the LVG results are consistent with the column densities 
derived assuming a constant rotational temperature. It is interesting to note that the H 2 densities derived 
from the molecular data are systematically smaller (a factor of ~ 3) than those derived from the dust 
observations (see Figure 17). 

The differences between the two determinations diminish in the low-density cores (a factor of ~ 2). These 
differences may partially be due to our assumptions concerning the dust millimeter opacity and temperature 
(see Galli, Walmsley, & Gongalves 2002 for a discussion of the latter and Bianchi et al. 2003 for a discussion 
of the former). Another possibility is that depletion of N 2 (and hence of N 2 H + and N 2 D + ) may occur in 
regions with density higher than 5 x 10 5 cm~ 3 (see also Caselli ct al. 1999; Bcrgin et al. 2001; Caselli et al. 
2002b; Crapsi et al. 2004; Belloche & Andre 2004). In the latter case, one would expect differences between 
maps of N 2 H+ and dust continuum (one possible example could be L183, see Section A). We note also that 
there is a hint in our data that densities derived from N 2 D + are systematically higher than those from N 2 H + 

Table 11 

Global and local velocity gradients from N 2 H+(l-0) maps. 



core 


G 


0" 


<Gi > b 




km s -1 pc -1 


O 


km s _1 pc -1 


L183 


1.4±0.1 


-49±1 


1.5±0.7 


OphD 


1.6±0.1 


-119±1 


1.6±1.0 


L694-2 


0.7±0.1 


-18±1 


0.8±0.3 


L429 


0.4±0.1 


-33±2 


1.9±1.1 


L492 


0.2±0.1 


83±2 


1.5±0.8 



a Direction of increasing velocity, measured East of 
North. 

b Mean values of the magnitude of local gradients 
and corresponding standard error. 
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Fig. 17. — H2 central density derived from dust observation compared with the determination from 
N2H+/N2D+ observations. The density inferred from the dust emission is always higher than the one 
derived from the molecular data, although the difference seems to diminish in the lowest density cores. 
Moreover, the estimates done starting from N2D + data furnish moderately higher H2 density, consistently 
with the idea that N2D+ is a better tracer of the dust emission. Data points from N2D" 1 " were slightly offset 
in the x-axis in order to show them better. 



and this again may be consistent with a hole caused by N 2 freeze-out. However, we stress that these trends 
need confirmation and that measurements with higher sensitivity and angular resolution are needed. 
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Table 12 

h 2 volume density and n 2 h + -n 2 d + column densities derived in the lvg approximation. 







N 2 H+ 




N 2 D+ 


core 


n(H 2 ) 


A r (N 2 H+) 


n(H 2 ) 


7V(N 2 D+) 
3 10 12 cm~ 2 




10 5 cm~ 


3 10 13 cm- 2 


10 5 cm- 


L1521F 


1.3 


3.4 


4.5 


1.6 


TMC2 


1.7 


1.0 


1.5 


1.5 


L1544 


1.5 


2.8 


2.8 


4.6 


L183 


1.0 


2.0 


3.5 


2.8 


OphD 


3.2 


0.6 


4.0 


2.6 


B68 


0.6 


1.2 






L429 


1.6 


3.1 


4.3 


4.9 


L694-2 


1.1 


2.1 


2.2 


3.3 



Note. — Errors in the LVG determinations are estimated 
around 30% from comparisons with CTEX determiantions. 
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